Table 1. Log of the RXTE and OSSE observations. The ASCA observation, near-simultaneous with the RXTE observation 3, is described in Section 2. 
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No. 


Date 


Time [UT] 


Exposure [s] 


Count rate [s ^] 


Time [UT] 


Exposure [s] 


Count rate [s 


Time [UT] 


Exposure^ [s] 


Count rate [s ^] 




(1996) 






(4-50 keV) 






(15-200 keV) 






(50-1000 keV) 


1 


May 22 


19:05:21-19:47:45 


2176 


7437±34 


17:58:41-18:17:53 


815 


253.7±1.9 








2 


May 23 


14:28:17-15:05:05 


1792 


6929±29 


14:28:33-14:53:05 


1136 


319.7±1.5 








3 


May 30 


07:57:53-08:35:29 


1872 


4796±22 


07:57:37-08:18:25 


850 


163.3±1.6 








4 


June 17 


01:45:37-01:54:25 
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3357±17 


01:45:21-01:54:57 


1062 


77.8±2.2 


01:25:06-02:08:47 
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8.25±0.61 
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June 17 


04:57:53-05:16:17 
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5113±25 


05:00:01-05:17:05 
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117.4±1.8 


04:30:47-05:16:40 


4598 


7.52±0.63 
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June 17 
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07:47:23-09:12:52 
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11.85±0.77 
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June 18 


03:23:13-03:34:25 
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2893±15 
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65.6±2.0 
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2670 


5.79±0.81 
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June 18 


06:47:45-07:00:17 
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54.0±1.9 


06:04:56-07:30:08 


3506 
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9 
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09:59:45-10:25:21 
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6011±29 


09:59:45-10:25:53 
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140.9±1.6 


10:00:15-10:41:46 


4235 


8.59±0.64 


10 


Juno 14-25 














15:55:58-13:13:20'' 


670501 


9.833±0.098 



''The sum of two individual detectors live times. 
''The start and end time on June 14 and 25, respectively. 



Table 2. The fit results of the non-thermal model to the data. Values in parenthesis denote parameters fixed during the fits. Other important fixed parameters are: fa = 10, 7, 
1.3, i = 45°, ylpe = 1.0. See Section 3.8 for discussion of the fitted values of A/x. 
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ABSTRACT 

We present X-ray/7-ray spectra of Cyg X-1 observed during the transition from the 
hard to the soft state and in the soft state by ASCA, RXTE and CGi?0/OSSE in 
1996 May and June. The spectra consist of a dominant soft component below ~ 2 
keV and a power-law-hke continuum extending to at least ^ 800 keV. We interpret 
them as emission from an optically-thick, cold accretion disc and from an optically- 
thin, non-thermal corona above the disc. A fraction / ^ 0.6 of total available power 
is dissipated in the corona. 

We model the soft component by multi-colour blackbody disc emission taking 
into account the torque-free inner-boundary condition. If the disc extends down to the 
minimum stable orbit, the ASCA/ RXTE daia. yield the most probable black hole mass 
of Mx ~ IOMq and an accretion rate, M « 0.5Le/c^, locating Cyg X-1 in the soft 
state in the upper part of the stable, gas-pressure dominated, accretion-disc solution 
branch. 

The spectrum of the corona is well modelled by repeated Compton scattering 
of seed photons from the disc off electrons with a hybrid, thermal/non-thermal dis- 
tribution. The electron distribution can be characterized by a Maxwellian with an 
equilibrium temperature of kT^ ^ 30-50 keV and a Thomson optical depth of r ~ 0.3 
and a quasi-power-law tail. The compactness of the corona is 2 ^ ^ 7, and a 
presence of a significant population of electron-positron pairs is ruled out. 

We find strong signatures of Compton reflection from a cold and ionized medium, 
presumably an accretion disc, with an apparent reflector solid angle, Q,/2tt ~ 0.5-0.7. 
The reflected continuum is accompanied by a broad iron Ka line. 

Key words: accretion, accretion discs - radiation mechanisms: non-thermal - stars: 
individual (Cygnus X-1) - gamma-rays: observations - gamma-rays: theory - X-rays: 
stars 



1 INTRODUCTION 

Cyg X-1 is the best-studied Galactic black-hole candidate. 
Its X-ray source was discovered in a 1964 June rocket flight 
(Bowyer et al. 1965). Its optical companion is HDE 226868 
(Webster & Murdin 1972; Bolton 1972), an OB supergiant 
(Walborn 1973). HDE 226868 belongs to the NGC 6871/Cyg 
0B3 association (Massey, Johnson & DeGioia-Eastwood 
1995). The distance to the association has been measured 



by Massey et al. (1995) to be D ~ 2.1 ± 0.1 kpc, whereas 
Malysheva (1997) obtained D ~ 1.8 kpc using a different 
photometric system. Hereafter, we adopt D — 2 kpc. We 
note that a commonly-used value of 2.5 kpc was obtained 
from the lower limit set by the extinction observed towards 
Cyg X-1 compared to the extinction-distance dependence 
for stars with directions close to that of Cyg X-1 (Bregman 
et al. 1973; Margon, Bowyer & Stone 1973). However, this 
method cannot be used for OB supergiants, in which case 



© 1999 RAS 



2 M. Gierlmski et al. 



a significant part of tlie extinction is commonly focal (e.g. 
Massey et af. f995). 

Tfie spectroscopic mass function is 



/m = 



Mx sin'^ i 



(0.252 ±0.010)Mo 



(1) 



(f + Mo/Mx)2 

(Gies & Bofton 1982), wliere i is ttie incfination of the orbitai 
axis to the observer, AIq and Mx are the masses of the 
supergiant and the compact object, respectiveiy. The mass 
ratio, q = Mq /Mx, has been constrained to be within 1.5 ^ 
q < 2.3 (Ninkov, Walker & Yang 1987; Gies & Bohon 1986). 
Thus, Mx is strongly anticorrelated with i, 



1.5M0 <, Mxsin I ^ 2.9Mg 



(2) 



The companion mass, Mq, is probably in a range between 
~ 15Mo (Herrero et al. 1995) and ~ SOM© (Gies & Bohon 
1986). This mass range together with the above range of q 
lead to Mx between ~ 6.5M0 and ~ 2OM0. Combined with 
the mass function, these constraints imply 25° ^ i ^ 50°. 
This is within limits from various considerations (e.g. Dolan 
& Tapia 1989; Ninkov et al. 1987; Davis & Hartman 1983) 
yielding 25° < i < 67° (where the upper limit is due to 
the lack of eclipses), and a recent Ha Doppler tomography 
(Sowers et al. 1998), giving i < 55°. 

The X-ray/7-ray (hereafter X7) emission of Cyg X- 
1 undergoes transitions between a hard (also called 'low') 
spectral state and a soft ('high') one. Most of the time the 
source is in the hard state, where the X7 spectrum can be 
roughly described by a hard power law with a photon index 
of r ~ 1.6-1.8 and a high-energy cut-ofl above ~ 100 keV, 
and a Compton reflection component including an Fe Ka 
line (e.g. Gierlihski et al. 1997a, hereafter G97). A soft ex- 
cess observed below 3 keV includes a blackbody with tem- 
perature fcTbb ~ 0.15 keV (Balucihska-Church et al. 1995; 
Ebisawa et al. 1996) . A likely model for the hard state of Cyg 
X-1 embodies Comptonization of blackbody photons from 
an optically-thick accretion disc in a hot, thermal, optically- 
thin plasma. 

In the soft state, the X7 spectrum is dominated by a 
soft, blackbody-like component peaking at 1 keV followed 
by a tail, which can be roughly described as a power law with 
r ~ 2.5 (Dolan et al. 1977; Ogawara et al. 1982), extending 
to at least several hundred keV (Phlips et al. 1996) . The soft 
component can be modelled as thermal emission from an 
optically-thick accretion disc (Mitsuda et al. 1984; Hanawa 
1989). The nature of the hard tail has remained unclear. 
The soft state of Cyg X-1 was observed in 1970, 1975, 1980, 
1994 and 1996, each time lasting no more than a few months 
(Liang & Nolan 1984; Cui et al. 1997). In 1994, the soft state 
was observed only in hard X-rays and soft 7-rays by the 
OSSE and BATSE detectors aboard Compton Gamma Ray 
Observatory (Phlips et al. 1996; Paciesas et al. 1997; Ling 
et al. 1997), with no accompanying observations in soft X- 
rays. Figure ^ shows X7 spectra of Cyg X-1 in the hard and 
soft states of 1991 and 1996, respectively, as well as in an 
intermediate state of 1996 May. 

In this paper, we study RXTE, ASCA and OSSE spec- 
tral data from observations of Cyg X-1 during its soft state 
in 1996. We find that spectra can be well described by emis- 
sion of a non-thermal plasma in the vicinity of an accretion 
disc. We study then implications of this model for the size 
of the plasma and its location with respect to the disc, the 
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Figure 1. Spectral states of Cyg X-1. The hard-state spectrum 
was observed by Ginga and OSSE on 1991 June 6 (G97) and by 
COMPTEL 1991 May 30-June 8 (McConnell et al. 1994). The 
intermediate state was observed by RXTE on 1996 May 23 (data 
set 2 in Table Q). The soft state was observed by ASCA and RXTE 
on 1996 May 30 (data set 3) and by OSSE 1996 June 14-25 (data 
set 10). Here and in subsequent figures, the spectral data are re- 
binned for clarity. The curves correspond to the hybrid model 
described in Section 3.4 



presence of e^e pairs, the inner radius of the accretion disc, 
the disc stability, and the mass of the central black hole. 



2 THE DATA 

The soft state of 1996 lasted approximately from May 16 
to August 11 (Zhang et al. 1997c). The observation log of 
RXTE and OSSE observations is given in Table |l[ 

The RXTE data, from the HEASARC public archive, 
were obtained during observations with proposal numbers 
10412 and 10512. After preliminary selection, we analysed 
the data of 1996 May 22, 23, 30, and June 17, 18. During 
the reduction, we rejected data affected by the South At- 
lantic Anomaly and those for which elevation above Earth's 
horizon was less then 10° . We then constructed PCA pulse- 
height spectra summing up counts from all PCA layers and 
detectors. We obtained HEXTE spectra separately for clus- 
ter A and B. We added a 2 per cent systematic error to 
each PCA channel to represent residual uncertainties of cal- 
ibration. Since dead-time effects of HEXTE clusters are not 
yet well established, we allowed free normalization of the 
HEXTE spectra during spectral fits. We used the PCA re- 
sponse matrices v3.0 as included in the ftools v4.1 released 
on 1998 May 8, and the HEXTE response matrices of 1997 
March 20. 

OSSE observed Cyg X-1 on 1996 June 14-25 (viewing 
period 522). The lightcurve of this observation is shown in 
Figure ^ The flux is moderately variable, within a factor of 
~ 2. The OSSE observation overlaps with 6 RXTE observa- 
tions on June 17 and 18 (data sets 4-9). We have thus ex- 
tracted 6 OSSE data sets (in the 50-1000 keV range) nearly 
simultaneous with the RXTE observations. However, in or- 
der to get better statistics, we have increased each OSSE 
data interval by 30 minutes on each side of the correspond- 
ing RXTE interval, as shown in Table §. The OSSE data 
include energy-dependent systematic errors (estimated us- 
ing both in-orbit and pre-launch calibration data), which 
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Figure 2. The OSSE light-curve in the 50-150 keV band. Each 
data point corresponds to one orbit. The boxes in the lower part 
of the plot indicate the times of six RXTE observations (Nos. 4—9 
in Table |]). 
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range from ~ 3 per cent at 50 keV to ~ 0.3 per cent at 
^ 150 keV. We have also obtained the OSSE spectrum from 
the entire observation (data set 10). The exceUent statistics 
of this spectrum allows a study of the average spectral prop- 
erties in soft 7-rays, e.g. possible presence of an annihilation 
feature or a high-energy cut-off. 

ASCA observed Cyg X-1 from 1996 May 30 5''30'"' to 
May 31 3''20™ UT (Dotani et al. 1997). The CIS observa- 
tion was made in the standard PH mode. The SIS data suf- 
fer from heavy photon pile-up and thus are not usable. The 
ASCA and RXTE observations on May 30 overlap only on 
a short period of time. Therefore, in order to improve the 
statistics, we used some ASCA data before and after the 
RXTE observation of the data set 3, which resuhs in 1488 s 
of the net exposure time (which corresponds to about 4000 
s of on-source data). Although the resulting ASCA/ RXTE 
data are not strictly simultaneous, the constancy of the 
ASCA light-curve, shown in Figure pi strongly suggests con- 
stancy of the ASCA spectrum. We used the GIS response 
matrix v4.0. 

Near-simultaneous ASCA/ RXTE observations of this 
bright source provide us with an opportunity to study the 
cross-calibration of the instruments. Since the response ma- 
trices of X-ray detectors are not diagonal, we cannot com- 
pare detector counts directly. We can only choose a model 
spectrum, fit it to the data, and then compare the resulting 
spectra. For this, we use our best model (consisting of a disc 
bla ckbody and a non-thermal component) found in Section 



3^ below. We first fit the 0.7-10 keV ASCA data simul- 
taneously with the 2-50 keV PCA and, in order to better 
establish the high-energy continuum, with the 15-200 keV 
HEXTE data. Figure ^ shows the data-to-model ratio for 
the best fit. We see that there are strong positive residuals 
at the high and low-energy ends of the GIS and PCA spec- 
tra, respectively, indicating that the PCA spectrum is softer 
than the GIS one in the range of the overlap. In particular, 
in the lowest-energy range of ~ 2.2-2.6 keV, the PCA chan- 
nel flux exceeds the corresponding GIS flux by around 30 
per cent, which corresponds to about 12a. This discrepancy 
is probably due to uncertainties in the PCA response ma- 



Figure 3. The 0.7-10 keV ASCA/GIS and 4-10 keV RXTE/PCA 
light-curves (time is measured from the beginning of 1996 May 
30) of the nearly simultaneous observations used here (data set 3 
in Table |). 



trix, especially at the lowest energies. A better agreement 
between the GIS and PCA data can be obtained if only the 
PCA channels at > 4 keV are used, as shown in Figure 
Then the PCA-to-GIS flux ratios at 4-4.3 keV and 9-10 keV 
become 0.96 and 1.08, respectively and we use the 4-50 keV 
PCA data hereafter. The above uncertainty in determin- 
ing the observed X-ray spectrum will then result in some 
systematic uncertainty in the parameters of our theoretical 
models. 

After completion of the data fits presented in this pa- 
per, an updated version of the PCA response matrices was 
released as a part of ftools v4.2 on 1998 December 10. We 
have then repeated the joint ^5Cj4/PCA fit above with the 
new PCA response matrix, and found that the discrepancies 
reported above have remained virtually unchanged. 



3 RESULTS 

3.1 Models of absorption, disc emission and 
reflection 

For spectral fits, we use XSPEC vlO (Arnaud 1996). The con- 
fidence range of each model parameter is given for a 90 per 
cent confidence interval, i.e., Ax^ = 2.7 (e.g. Press et al. 
1992). The spectra are absorbed by the interstellar medium 
with a column density A^h , for which we use the opacities of 
Morrison & McCammon (1983) and the abundances of An- 
ders & Ebihara (1982). Since only the ASCA data constrain 
A^'h (see Table H below), we assume A'^h = 5 x 10^^ cm~^ in 
all fits not including the ASCA data. 

We fit the observed spectra by a sum of 2 main primary 
continua. The first one is due to blackbody emission of an ac- 
cretion disc. We model it using the pseudo-Newtonian (here- 
after PN) potential and taking into account the torque-free 
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Figure 4. The best fit of tfie non-tliermal model to the near- 
simultaneous ASCA and RXTE data of May 30 (observation 3). 
The upper and lower panels show the data-to-model ratios for 
ASCA/GIS and RXTE/'PCK, respectively, (a) 0.7-20 keV resid- 
uals of the joint fit to the 0.7-10 keV ASCA data and the 2-50 
keV PCA data, (b) The residuals to the fit with the PCA data 
restricted to the 4-50 keV range. Note a resulting significant im- 
provement of the joint fit. 



in a strong gravitational potential, Doppler and gravita- 
tional shifts become important. We approximate these ef- 
fects convolving both the reflected component and the Fe 
Ka line with the Schwarzschild disc line profile of Fabian et 
al. (1989). For simplicity, the emissivity of the reflection and 
of the line as a function of a disc radius is assumed to vary 
as . We expect s ~ 0, 2, and 3 from different parts of the 
disc, depending on the geometry (see Fabian et al. 1989). 
We assume s = 2 unless stated otherwise. In calculating the 
relativistic distortion, we consider a range of radii from rin, 
which equals that used for the disc blackbody emission, to 
rout = 1000. 

The reflected medium is photo-ionized, and the ioniza- 
tion is described by a parameter ^. To allow direct com- 
parison with other results, we follow Done et al. (1992) in 
the definition of the ionization parameter as ^ = 47rFion/n, 
where Fion is the 5 eV-20 keV irradiating flux in a power- 
law spectrum and n is the density of the reflector. (We note 
that in the soft state of Cyg X-1, most of the ionizing pho- 
tons have energies J; 1 keV, and a power-law shape of the 
irradiating continuum is only a rough approximation.) For 
simplicity, we assume that f is uniform in the medium. Since 
the ionization state of the disc close to the black hole will 
vary significantly with radius, the value of ^ found from the 
fits corresponds to an average ionization of matter at a range 
of radii from which most of reflection originates. The abun- 
dances in the reflector are of Anders & Ebihara (1982) ex- 
cept the relative Fe abundance, Apc, which is allowed to be 
free in some models. The ion edge energies and opacities 
are from Reilman & Manson (1979) except that now the Fe 
K-edge energies taken from Kaastra & Mewe (1993). 

In our fits, we allow for free relative normalization of 
data from different instruments, except for the simultane- 
ous PCA/OSSE data. In that case, we found their relative 
normalization is fully consistent with unity. 



inner boundary condition, see Appendix^. The disc extends 
from an inner radius, i?in, corresponding to rin = i?in/i?g 
(where 7?g = GM/c?) to infinity, and its temperature distri- 
bution is parametrized by the maximum temperature, Tmax 
reached by the disc. We assume nn = 6 (i.e., at the mini- 
mum stable orbit in the Schwarzschild metric) unless stated 
otherwise. We assume D = 2 kpc and, in most fits, the disc 
inclination of i = 45° (which is around the middle of the 
allowed range of i, see Section |l|) and the colour-to-effective 
temperature ratio, /coi = Tcoi/rcff ~ 1-7 (Shimura & Taka- 
hara 1995). 

The second primary continuum corresponds to the ex- 
cess emission observed above the high-energy part of the 
disc emissio n. For t hat continuum, we use models described 
This continuum is emitted outside the 



in Sections 3.2 



3.4 



disc, and thus it is Compton-refiected from the disc surface. 

For modelling Compton reflection, we use viewing- 
angle-dependent Green's functions of Magdziarz & Zdziarski 
(1995), which assume an isotropic point source (or, equiva- 
lently an optically thin corona) above a slab. However, we 
treat fl/2n as a free parameter, where il is an effective (i.e., 
corresponding to the observed strength of reflection) solid 
angle subtended by the reflector. The reflection is accom- 
panied by an Fe fluorescence Ka line centred at an energy, 
Eiinc- When the reflection comes from a fast rotating disc 



3.2 Phenomenological description 

As illustrated on Figure 0, the X7 spectrum of Cyg X-1 in 
the soft state consists of a blackbody-like soft component 
dominant below ~ 2 keV and a power-law tail extending 
to at least several hundred keV. In order to allow a com- 
parison with previous results (Dotani et al. 1997, Zhang et 
al. 1997b) as well as to establish a phenomenological de- 
scription of the spectral data, we start our analysis with a 
simple model consisting of a PN disc and a power law with 
its corresponding Compton-reflection continuum and an Fe 
line computed taking into account the relativistic smearing. 
First, we fit the simultaneous RXTE/OSSF, data of June 
17 and 18 (data sets 4-9). The photon spectral index, F, 
varies between about 2.3 and 2.5, which confirms the re- 
sults of simultaneous ASM and BATSE monitoring of Cyg 
X-1 (Zhang et al. 1997b). The refiector solid angle, ^l/2n, is 
~ 0.6-0.7, and xl ~ 500/475. 

Next, we test the RXTE/OSSF, data for a presence of a 
high-energy cut-off. We multiply the power law by an expo- 
nential factor with an e- folding energy, E{. Since the simulta- 
neous RXTE/OSSF, observations constrain E{ only weakly, 
we test the properties of the high-energy tail using the OSSE 
data integrated over the whole viewing period 522 (data set 
10), which spectrum is detected up to ~ 800 keV (Figure 
When we fit together the RXTE data No. 5 and the inte- 
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Figure 5. OSSE spectra of Cyg X-1 in the soft state of 1994 
February (open boxes) and 1996 June (filled boxes). The former 
data come from the OSSE viewing period 318, 1994 February 1- 
8 (Phlips et al. 1996). The latter data are described in Table ^ 
observation 10. The two data sets are fitted by a power law with 
r = 2.72^0 and P = 2.57lo°3, respectively (solid lines). 



grated OSSE spectrum, the best-fit e-folding energy is very 
high, Ei = 940+^°° keV. 

On the other hand, a power law hard continuum fails 
to reproduce the low-energy spectrum of the simultane- 
ous ASCA/ RXTE observation, which was pointed out by 
Gierlihski et al. (1997b) and C98. A fit of a model consist- 
ing of a blackbody disc, a power law and the corresponding 
reflection continuum with an Fe Kq line provides a very 
poor fit to this data set, with xl ~ 1100/576. A much bet- 
ter fit, with xl = 608/574, can be obtained with a broken 
power-law continuum instead of a single power law. We find 
the power- law indices, Ti = 2.67to°ol, ^2 = 2.19 ± 0.02, the 
break energy, Streak = (10.9 ± 0.5) keV, D./2n = OASto.ol, 
Nu = (5.8 ± 0.1) X 10^^ cm~^ fcT^ax = 372t^ eV, Mx = 
9.4lg:^Mo, C = (2000 ± 600) erg cm s'^, Eune = 6.38l2:l5 
keV, and the line equivalent width, EW « 110 eV. These 
results are relatively similar to those of Dotani et al. (1997) 
and C98. 

Furthermore, we stress that if the tail beyond the disc 
blackbody is due to Comptonization of the disc blackbody 
photons (as it is most likely the case in Cyg X-1), the tail 
spectrum would have a sharp low-energy cut-off at 2 keV. 
A proper description of the tail is then essential for deter- 
mining the disc parameters. Using of a continuum without 
an intrinsic low-energy cut-off (e.g. a power law) can result 
in an erroneous determination of the disc properties and A'^e . 



3.3 Thermal plasma model 

As shown in the previous section, the high-energy power 
law in the soft state extends to very high energies {Et ^ 600 



keV). This contrasts the hard state, where the power law 
is cut off above ~ 100 keV. If this spectrum were emitted 
by Comptonization of soft photons in a thermal plasma, 
the electrons would have to have a temperature of at least 
several hundreds keV. Then, in order to account for the ob- 
served power-law slope, the optical depth of the electrons 
would have to be very low, ^ 0.01. Subsequent scattering 
events in a plasma with such very low optical depth would 
form distinct humps in the emerging spectrum, which are 
clearly not seen in the data. We confirm this using a thermal- 
Comptonization model of Coppi (1992), assuming that seed 
photons for Compton up-scattering are those emitted by the 
disc. This model (together with a disc, reflection and a line) 
fitted to the RXTE/OSSE spectra (data sets 4-9) yields un- 
acceptable values of ~ 1000/476, much worse than that 
obtained in section with a power-law incident continuum. 

Furthermore, we find that the thermal- Comptonization 
model provides a very poor fit to the ASCA/ RXTE data, as 
it is unable to reproduce the broken power-law shape of the 
tail. This result agrees with that of Cui et al. (1998, hereafter 
C98) who could fit those data with thermal Comptonization 
only assuming that the temperature of seed photons is much 
higher than that of the observed disc photons. Thus, we 
rule out the thermal Comptonization as a model for the 
continuum emission in the soft state of Cyg X-1. 



3.4 Hybrid thermal/non-thermal plasma model 

Since the high-energy tail in the spectrum cannot be due 
to emission from a thermal plasma, we fit the data using a 
non-thermal model of Coppi (1992) with some modifications 
(see also Coppi & Blandford 1992; Zdziarski, Coppi & Lamb 
1990; Zdziarski, Lightman & Maciolek-Niedzwiecki 1993 and 
a recent review by Coppi 1999). This model embodies Comp- 
ton scattering, pair production and annihilation, pe* 
and e^e* thermal and non-thermal bremsstrahlung, and en- 
ergy exchange between thermal and non-thermal parts of the 
distribution via Coulomb scattering. The seed photons 
for Compton scattering are those emitted by the PN disc. 
Selected e~ are accelerated to suprathermal energies and 
the thermal part of the distribution can be additionally 
heated (apart from Compton heating/cooling and Coulomb 
heating by the non-thermal e*). [A similar, thermal/non- 
thermal, model was used, e.g., by Holman & Benka (1992) 
for solar flares.] 

This model assumes a spherical plasma cloud with 
isotropic and homogeneous distributions of photons and e*, 
and soft seed photons produced uniformly within the plasma 
(Coppi 1992). The properties of the plasma depend on its 
compactness, £ = Ca^ / (TZnicC?), where yC is a power of 
the source, TZ is the radius of the sphere and is the 
Thomson cross section. We specifically use a hard compact- 
ness, ^h, which corresponds to the power supplied to the 
electrons, and a soft compactness, ^s, corresponding to the 
power supplied in the form of soft seed photons. We as- 
sume TZ = 10* cm, but this value is important only for 
bremsstrahlung emission (which is negligibl e in our case), 
and for the Coulomb logarithm (see Section 4.2 below). 

This spherical model is certainly an idealization of the 
hot plasma source in Cyg X-1, and we need to relate the 
model parameters to those corresponding to the actual X7 
source in Cyg X-1. Based on the strength of the blackbody 
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Figure 6. A sketch of the disc-corona geometry. The emission of 
the spherical plasma model is identified with that of the spherical 
region shown here. Power dissipated in the disc and corona arc 
Pd and Pc, respectively. Lg, Lh, ir a re s oft, hard and reflected 
luminosities, respectively. See Section and Appendix B. 



and reflection components (Section 3.2), our preferred ge- 
ometry is that of a hot corona above the surface of an ac- 
cretion disc. Similar to a method of Bjornsson & Svensson 
(1991), we consider a local spherical volume with a radius 
TZ = H/2 inside the corona, where H is the corona scale- 
height, as shown in Figure H. In some models, we will also 
allow the corona to cover the disc at ii < 7?h partially with 
the covering fraction, g (< 1). 

We then consider the relation of the compactnesses, 
and £s, obtained from fitting the spherical model, to the 
parameters of the corona: the luminosity in scattered pho- 
tons, Lh, the luminosity in unscattered blackbody photons, 
Ls, the characteristic radius of the corona, i?h, and H. This 
will allow us to constrain the last two quantities by the ob- 
served fiuxes and fits of the model. For that purpose only, 
we assume that the corona radiates isotropically and ne- 
glect corrections due to the finite Thomson optical depth of 
the corona, the finite disc albedo, a, and the viewing angle 
of i 7^ 60° (in which case the anisotropy of the blackbody 
fiux emitted by the disc plays a role). However, we do take 
into account those corrections later in considering the en- 



3.9 



ergy balance between the disc and the corona (Section 
and Appendix B). 

We point out here one significant difference between the 
spherical and disc-corona geometries. Namely, the sphere 
radiates away all of the power supplied to the electrons as 
well as that supplied as seed photons. On the other hand, the 
isotropic corona luminosity, Lh, corresponds to only half of 
the total coronal power, 2Lh (while the other half is radiated 
towards the disc and absorbed), whereas the disc luminosity, 
Ls, corresponds to all of the power in soft photons emitted 
by the disc. In our model fitting, we fit observed fiuxes in soft 
and hard photons. Therefore, we have to define in terms 
of the hard power radiated away (rather than the dissipated) 
by electrons in the spherical volume. This yields 



We note that the numerical factor in this equation is not pre- 
cise due to the approximations described above, which will 
introduce some systematic uncertainty in our determination 
of plasma parameters. 

In models in this Section, we assume g = 1, which cor- 
responds to all the blackbody photons passing through the 
plasma. Then, the soft compactness. Is, is to be defined in 
the same manner as in equation (^) except that Lh is re- 
placed by Ls. If g < 1, we will distinguish between the soft 
compactness in photons incident on the hot plasma and that 
in soft photons freely escaping (Section 3.9). 

The compactnesses corresponding to the electron accel- 
eration and to the additional heating of the thermal part 
of the distribution are denoted as ^nth and £th, respec- 
tively, and £h = Inth + £th- In our fits, we use the ratio of 
^nth/^h as a parameter. The rate at which non-thermal elec- 
trons appear in the source is assumed to be a power law, 
iVinj(7) c>c 7~'"'"j, between the Lorentz factors of 7min and 
7max. We hereafter assume jrnin = 1-3 and 7max = 1000. A 
different choice of 7max has a negligible effect on our results 
as long as 7max 3> 7min siuce most of the power in the case 
of Finj > 2 (see Table 2) is injected around 7min. The effect 
of t he choice of 7min on our results is discussed in Section 
l2[ 

The injected electrons and pairs produced in photon- 
photon collisions lose energy by inverse- Compton scattering 
(at a rate of 7compton oc 7^^), so the equilibrium distri- 
bution is softer then that of injected electrons. This can be 
seen from the steady-state continuity equation. 



Mnj(7), 



(4) 



which solution leads to the non-thermal component of e 
distribution, Nc{'y) oc 'y~('"inj + ^). Therefore, the single scat- 
tering of seed photons off these forms a power-law spec- 
trum with photon index F « Finj/2 -I- 1. In steady state, the 
non-thermal and thermal have optical depths, Tnth (typ- 
ically <C 1) and Tth, respectively, following from balance of 
acceleration, pair production, and energy losses. The total 
optical depth is denoted by r. 

The optical depth of the thermal component, Tth (typi- 
cally 3> Tnth), is a sum of the optical depth to scattering on 
e'^e~ pairs and on e~ coming from ionization of atoms, ri. 
The latter one is a free parameter of the model, but nh as 
well as the equilibrium temperature of the thermal compo- 
nent, Tc, are computed self-consistently. 

The spectrum of the model is then used to compute the 
corresponding continuum fro m Co mpton reffection including 
relativistic smearing (Section 3.1). We find that the plasma 



in the source is so energetic that a single Compton scat- 
tering effectively removes a Compton-reffected photon from 
the reflection spectrum. Then, fl/2n obtained from the fit 
corresponds to the unscattered part of Compton reflection. 

A sequence of model spectra for varying compactness is 
presented on Figure ^ For parameters consistent with the 
observations (as found in Section 3.5 below) the spectrum is 
basically a power law without a break up to ~ 10* keV. The 
power law is formed by a single inverse- Compton scatter- 
ing of disc seed photons off the non-thermal with optical 
thickness of Tnth ~ 10~^. At low energies, we see the spec- 
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Figure 7. Model spectra from non-thermal plasmas above the 
surface of an accretion disc calculated for varying soft compact- 
ness, 4- Prom bottom to top, 4 = 0.01, 0.1, 1, 10, 100 and 1000. 
Other parameters are £h/^s = 0.3, ^nth/^h = 0.9, Ti = 0.3, Finj = 



3, 7rr 



1.3, fcTn: 



: 0.4 keV. At low and high compactness, we 



see distortio ns o f the hard spectrum due to Coulomb interactions 
(see Sec tion |4.2| ) and e+e~ pair production and annihilation (see 
Section |3.6| ), respectively. Compton reflection and seed photons 
escaping from the plasma without a scattering are not included 
in the spectra shown. The normalization of each spectrum is ar- 
bitrary. 



trum of disc photons repeatedly scattered off the thermal 
e*, as well as singly-scattered off the nonthermal e*. The 
former scattering dominates, forming a tail beyond the disc 
blackbody energies, crossing the non-thermal power law at 
~ 10 keV. 



3.5 Spectral fits with the hybrid model 

In fitting the data, we have found we cannot uniquely deter- 
mine the overall level of the compactness, which, in our fits, 
is given by the value of Is- The joint constraints from both 
the ASCA/RXTE and RXTE/OSSY. data yield 5^4^ 20. 
Since the ratio of 4/4 = O.Set^ Jlla 

is well established (see 
below), this corresponds to constraints on the hard com- 
pactness, 2 ^ 7. The upper limit is due to the onset 
of efficient e^e~ pair production at ^ 7 (see Figure |^, 
whereas the data do not show any signatures of the pres- 
ence of pairs (see Section ^^). The lower limit is due to 
Coulomb cooling of non-thermal electrons becoming impor- 
tant at low compactness (see Figure |^ . At small this pro- 
cess thermalizes efficiently the n on-t hermal electrons below 
some Lorentz factor (see Section ^!2| ) , whereas the data are 
consistent with a non-thermal distribution close to a power 
law. We hereafter assume Is = 10 in most of fits below, and 
use lii/ls rather than as a free parameter. 

We start spectral fitting with the June 17-18 
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Figure 8. The simultaneous RXTE and OSSE observation of 
Cyg X-1 in the soft state on 1996 June 17 (observation No. 5) 
together with a model with t\ = 0.25 (solid curve; fit parameters 
given in Table The model is decomposed into the hot-plasma 
emission (short dashes), Compton reflection including the Fe Ka 
line (dots), and the disc blackbody emission (long dashes). 



7?Xr£/0SSE data (observations 4-9). We find that the data 
do not allow us to determine ti. There are several shallow 
local minima of in a wide range of Ti between ~ 0.2 
and ~ 3. However, models with different values of Ti dif- 
fer significantly from each other below 4 keV, where we 
have no data. We find a strong correlation between Ti and 
Ih/ls, which latter parameter is also poorly determined by 
the data. Therefore, in order to determine Ti and (.ti/ls, we 
use the ASCA/RXTE observation, for which we obtain Ti 
= 0.25 



+0.05 
-0.04 



and 4/4 = 0.36: 



(see below). Since the 



spectral properties did not change significantly between 1996 
May 30 and June 17-18, we assume that Ti also remained 
similar, and we fix Ti = 0.25 in subsequent fits to the data 
sets 4-9. Figure ^ shows the data set 5 fitted with Ti = 0.25. 
We further note that a large value of Ti, e.g. 2.5, would be 
difficult to reconcile within the coronal model, which we find 
below the most likely, with the presence of strong Compton 
refiection in the spectrum. 

For the data sets 4-9, we find the best-fitting values of 
^nth/4 ~ 0.95-1 and corresponding 90 per cent confidence 
intervals covering the range of ~ 0.8-1. Thus, in order not to 
introduce unnecessary free parameters, we assume ^nth/^h = 
1 in our fits to the TJXTE/OSSE data (see Table 1]). 

Subsequently, we study the ASCA/RXTE data of May 
30. The 0.7-10 keV energy range of the GIS data allows us to 
determine the A^'h, and therefore we let it free in this fit. The 
low-energy coverage also allows us to better constrain the 
direct heating of thermal electrons, and thus we let £nth/^h to 
be free. The fit results are shown in Table ^ The best-fitting 
parameters correspond to Tth ~ 0.26, kTc « 46 keV, and 
Tilth ~ 10~^. We find now the relative non-thermal power 
of 4th/4 = 0.77ini 

which, although somewhat less than 
unity, still confirms our conclusion that the direct thermal 
heating in Cyg X-1 in the soft state is weak compared to 
non-thermal acceleration. The above range of £nth/^h is still 
consistent with the resuhs of June 17-18 iJXTE/OSSE fits 
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Figure 9. The simultaneous ASCA and RXTE observation of 
Cyg X-1 in the soft state on 1996 May 30 (observation No. 3). 
The short dashes, dots and long dashes show the Comptonization 
continuum, the component reflected from the cold disc and the 
disc emission, respectively. The solid curve shows the sum. The fit 
parameters are given in Table ^. The Comptonization continuum 
is dominated by thermal and non-thermal scattering below and 
above, respectively, ~ 15 keV, with the resulting spectral break 
at that energy. 



within the error bars, and thus we are not able to conclude 
whether the relative fraction of non-thermal power changed 
between May 30 and June 17-18 or not. The range of Is 
found above corresponds, at the fitted ^h/^s ~ O.SGIqqjj ^o 
the total compactness, I = Is + (-h, of 7 ^ £ ^ 27. 

We note here that the actual value of ^nth/^h depends 
somewhat on the assumed value of Is . This effect can be un- 
derstood based on results of Section below. Namely, the 
importance of Coulomb heating of the thermal plasma 
component by non-thermal electrons decreases with increas- 
ing compactness. Thus, more direct heating of the thermal 
is required at a higher compactness than at a lower one. 
Indeed, for the lower and upper limit of Is allowable by the 
data, 5 and 20, we find = 0-S8to°ol and 0.72t°;°^, 

respectively. 

On 1996 May 22-23 (data sets 1 and 2), Cyg X-1 was in 
an intermediate state. The hard X-ray fiux was significantly 
stronger than in soft-state observations of May 30 and June 
17-18, but weaker than in the hard state (see Figure |l|). 
Although Cyg X-1 was observed during this period only by 
RXTE, we find that the shape of the spectrum allows a 
determination of ri . Fit results are given in Table ^. Since the 
soft component is relatively weak (see Figure |l]), we can only 
constrain the seed photon temperature, fcTmax, to be ^ 0.4 
keV, and we thus assume fcTmax = 0.3 keV. We caution here 
that fcTmax is strongly anti-correlated with _Rin. For example, 
the data set 1 at fcTmax = 0.15 keV yields Rin about 6 times 
larger than that for fcTmax ~ 0.3 keV. We also obtain the 
compactness ratio of i^/is '~ 1, which indeed places these 
data between the hard and and soft states, for which £h/£s 
~ 10 and 0.3-0.4, respectively. 

Summarizing our results so far, we have found a satis- 
factory physical description of the data by a model with the 
main component due to Compton scattering of soft photons 



from an optically-thick accretion disc in an optically thin, 
hot, plasma. We find the plasma cannot be purely thermal. 
Rather, the electrons in the plasma have a hybrid distribu- 
tion with a Maxwellian component at a fcTo ~ 40-50 keV 
and a quasi power-law, non-thermal, tail. 



3.6 



pairs 



The optical depth of the thermal part of the e distribution, 
Tth, includes contributions from both ionization electrons 
(= Ti) and e+e" pairs. We find first that the ASCA/ RXTE 
observation (data set 3) cannot constrain the presence of 
pairs. Fits with Is ~ 40 and 10 with a difference of only 
2.2 yield rth/Ti 1 (a pair-dominated plasma) and 1.07 
(an e~ dominated plasma), respectively. This is because a 
presence of pairs manifests itself primarily by an annihilation 
feature around 511 keV (see Section 4.2), which energy range 
is not covered by RXTE. 

On the other hand, the RXTE/OSSB observations 
(da ta sets 4-9) do not constrain ri, as discussed in Section 



3.5 above, which again prevents us from constraining rth/ri. 
Therefore, we fit the ASCA/ RXTE data set 3 together with 
the average OSSE spectrum (data set 10). The presence of 
pairs is now excluded at a high statistical significance. At Is 
= 10, Tth/Ti « 1, and at Is = 40, which yields Tth/n « 1.4, 

becomes significantly worse, Ax^ ~ +16 at 626 d.o.f. The 
plasma becomes pair-dominated, rth/ri 3> 1, at i?s = 60, but 
then Ax^ = +20. The main contribution to the increasing 

comes from the OSSE channels around moC^ due to an 
annihilation feature produced by the model but not present 
in the data. Thus, pairs appear to be negligible. 

In our plasma model, we have assumed so far that the 
accelerated (or injected) particles are electrons rather than 
e^e~ pairs (and the acceleration rate is balanced by the 
rate at which electrons move from the non-thermal distri- 
bution to the thermal one). On the other hand, some physi- 
cal mechanisms may give rise to production of non-thermal 
pairs rather than electrons, as well as mostly pairs would be 
accelerated out of the thermal distribution if that is already 
pair-dominated. To investigate that case, we have consid- 
ered a plasma model with an injection of e^e~ pairs. How- 
ever, this model gives an extremely poor fit to the data, 

~ 1000/626. This is due to the luminosity from annihi- 
lation of the injected pairs, which equals 



-'annihilation 



Tinj - 2 
(Finj - l)7n 



(5) 



The annihilation luminosity is then large in our case with 
Finj ~ 2.5 (Table 2) and 7min constrained to be ~ 1 (see 
Section 4.2 below), which is in strong disagreement with the 
OSSE data showing no pair annihilation feature. Summariz- 
ing, the presence of a substantial number e^e~ pairs in the 
X7 source of Cyg X-1 in the soft state is ruled out. 



3.7 Compton reflection and Fe Ka line 

We find manifest signatures of reprocessing by a cold 
medium - an Fe K edge with a Ka line and a broad excess 
over the continuum at ~ 10-100 keV - in all our data sets. 
The most straightforward interpretation of those features is 
in terms of Compton refiection of the hot-plasma emission 
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by an ionized, optically-thick, accretion disc. The statisti- 
cal significance of the presence of reflection is very high. 
Fits to the data sets 3-9 without reflection are much worse, 
~ 2, than those including reflection (see Table ^), and 
the main cause of the bad fit are strong residuals around the 
Fe edge/line energies. The intermediate-state RXTE data 
sets, 1-2, are somewhat less sensitive to reflection. Still, fits 
without reflection are much worse [X^ = 654 and 586 at 425 
d.o.f.) than those in Table ^. 

Below, we use the simultaneous ASCA/ RXTE obser- 
vation (data set 3), which allows the most detailed study 
of Compton reflection and the associated Ka fluorescence. 
This is due to a broad energy range, ~ 0.7-100 keV, to- 
gether with good sensitivity near 7 keV. First, we confirm 
that the reflected continuum is accompanied by an iron Ka 
line, by means of the standard F-test (e^. Mises 1964) . We 
apply here our best model (see Table p[) with relativistic 



Compton reflection with a line (described in Section 3.1) 



and without it. Our null hypothesis is that a model with the 
line flts the data better than a model without it. A fit with- 
out the line yields Xi = 652 with v\ — 576 d.o.f., whereas 
a fit with the line yields xl = 618, V2 = 574 d.o.f. The 
chance function, ~ / lS.v)/{x^/v2) has a Snedecor's 
F distribution. The significance level, P{F > Fm), is the 
probability of committing an error of first kind, i.e. reject- 
ing the correct hypothesis. For the fits above, F^ = 15.6 
and P{F > Fm) = 2.5 x 10~^, which represents a very high 
statistical significance that an Fe Ka line is indeed present 
in the data. 

We then investigate the issue of the width of the line. 
For that purpose, we first compare the ASCA/ RXTE data 
to a model with static refiection (i.e., ignoring the relativistic 
effects) and a narrow Gaussian line. This yields a fit worse 
by Ax^ = +9 compared to the fit in Table ^. This means 
that indeed the data are much better described by the line 
and reflection being relativistically smeared than by static 
line and reflection. Still, the disc inner radius is only weakly 
constrained by relativistic smearing, to rin ^ 50 (in which fit 
we assumed the inner radius of the blackbody disc emission 
constant). Also, the value of the power-law index of refiec- 
tion emissivity (Section is only weakly constrained to 



2.1 



-1-0.5 
-0.7- 



In order to compare with results of C98, we 
also considered a broad Gaussian line in the static-reflection 
model, obtaining a rather large width of crpe = 1-18^q'22 keV 
(somewhat more than the values given by C98). Since uncer- 
tainties in the RXTE/PCA response matrix are not yet fully 
understood, we have checked our results by ignoring PGA 
data below 10 keV, in which case the Fe line parameters 
were constrained by ASCA data only. The results for the 
line energy and width turn out to be very similar to these 
obtained from the overlapping ASCA/ RXTE data. Thus, 
our conclusion is that the Fe Ka line in the soft state of Cyg 
X-1 is broad. 

The reflecting medium is ionized, with ^ = 350ti2o ^^S 
cm s~^. Though observations on May 22-23 (data sets 1-2) 
and June 17-18 (data sets 4-9) yield higher values of ^ (see 
Table ^ , those values are relatively weakly constrained and 
most of the data are consistent with the ASCA/ RXTE result 
within 90 per cent confidence. Note that this ionization, ob- 
tained with our physical incident continuum, is much weaker 
than that found with our phenomenological, broken-power 
law, continuum in Section 3.2, which shows the sensitivity 
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Figure 10. The dependence between the solid angle subtended 
by the reflector, Q/2it, and the assumed inclination angle, i. The 
error bars are given for 90 per cent confidence. 



of the fitted ionization parameter to the shape of the con- 
tinuum around 7 keV. 

The ionization state found here corresponds to a distri- 
bution of Fe ions peaking at Fe xxill^J,, and with Fe xxiv- 
XXVI constituting ~ 30 per cent of the ion population. Due 
to resonant trapping followed by the Auger effect (e.g. Ross 
& Fabian 1993), Ka photons following ionization of Fe xvil- 
XXIII ions are completely destroyed and yield no contribu- 
tion to the observed line. We thus would expect only a line 
from Fe xxiv-xxvi at ~ 6.7-7.0 keV. Instead, our model 
yields the rest-frame centre energy of i?iinc ~ 6.37 ± 0.14 
keV for the ASCA/ RXTE data and ~ 6.3-6.5 keV for other 
RXTE observations. This inconsistency is likely to be due to 
our simplified model of refiection, with the ionization state 
uniform throughout the refiecting medium, whereas the ion- 
ization would, in general, vary considerably with the disc 
radius. 

We also consider the Fe abundance, Afc, of the reflector. 
In all flts so far, we have assumed Apa = 1. We find that 
indeed Apa = 1 is consistent with the data. In particular, 
the ASCA/ RXTE data fitted with as a free parameter 
yield Afc = 1.0l°;t 

One parameter of Compton reflection which remains 
relatively uncertain is the inclination angle, i. As dis- 
cussed in Section |^, the most conservative range of i is 
25° < j < 67°. The resulting Q/2-k ~ 0.5-1.0, increas- 
ing with the increasing i, is shown in Figure |l^. At the 
Doppler-tomography limit of i = 55° of Sowers et al. (1998), 
n/2TT = 0.73 ± 0.11. Our fits favour large inclinations an- 
gles: = 626, 618, 612 (for 574 d.o.f.) at i = 25°, 45°, 67°, 
respectively, similarly to a recent result of Done & Zycki 
(1999) for the hard state of Cyg X-1. On the other hand, the 
plasma and disc properties, e.g., the values of £h/is, ^nth/^h, 
Ti, Tmax, are only weakly dependent on i. 



3.8 Inclination and the black-hole mass 

If we assume the disc extends to the minimum stable orbit in 
Schwarzschild metric, nn ~ 6, the ASCA/ RXTE data yield 
Mx = 13.3to gM0 at the assumed D = 2 kpc, f^oi = 1.7 and 
i = 45° . We note that this mass range is above that allowed 
by equation (^ at i = 45°. However, as found in Section 
a/A the model parameters other than Q/2tv are insensitive 
to the assumed value of i. T hus, A/x can be rescaled to other 
values of i using equation (AlO), which yields 
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11.2+^;°M 



(6) 



which holds for our data to an accuracy better than 2 per 
cent. This range of Mx overlaps with that of equation 
provided, 



0.13 < 



(1 -sin^ j)i/4 



< 0.28, 



(7) 



which is satisfied for 29° i ^ 39° , corresponding [equation 
(P)] to 14 ^ Mx/M0 ^11. These are the mass and inclina- 
tion ranges both consistent with fitting the X-ray data with 
a disc extending to the minimum stable orbit and satisfy- 
ing the constraints on the mass function and the mass ratio 
presented in Section |l|. 

However, this estimate neglects GR effects in the vicin- 
ity of a black hole, namely gravitational redshift and focus- 
ing. We can estimate these effects by introducing two fac- 
tors, /gr and (?GR, correcting the colour temperature and 
the integrated flux, respectively (Zhang, Cui & Chen 1997a; 
Cunningham 1975). Then, the correction factor to the mass 
for a non-rotating black hole, rin = 6, and i ~ 35° to the 
mass is QqI^^ fan ~ 0.83 (with a slow dependence on i). 
Thus, the above mass range becomes 



Mx f» (10± 1)M(7 



(8) 



which represents our best mass estimate including GR cor- 
rections. We note that this agrees very well with the esti- 
mates of Dotani et al. (1997) and C98 obtained using a GR 
disc model (after rescaling from D — 2.5 kpc used by them 
to 2 kpc). 

This estimate is sensitive to uncertainties in the dis- 
tance and the factor /coi, with Mx (x -D/coi- E.g., a 20 
per cent uncertainty in -D/coi would increase the range of 
i and Mx (after the GR correction) to 27° ^ i ^ 42°, 
14 ^ Mx/Mq ^ 7. 

Above, we assumed rin = 6. However, the data do not 
require that. When we allow rin to be a free parameter, 
the shape of the disc spectrum changes, and we find this 
change large enough to constrain rin by the ASCA/ RXTE 
data. Namely, the model fit to the data is improved by 



3.9 Energy balance 

If the hot plasma (either homogeneous or patchy) forms a 
static corona above the disc, energy balance can be used to 
obtain the fraction of the gravitational energy released in 
the corona, / (Svensson & Zdziarski 1994, hereafter SZ94). 
The general solution to this problem is given by equation 
(B8) in Appendix B. We assume here isotropic emission 
of the corona (corresponding to d = 1 in that equation). 
The reflection albedo, a, is computed to be a ~ 0.25 us- 
ing our ASCA/ RXTE model spectrum at i = 65°, at which 
the strength of Compton reflection is approximately equal 
to the reflection strength averaged over the 2n solid angle 
(Magdziarz & Zdziarski 1995). Since the plasma is optically 
thin, we assume the scattering probability, psc ^ Ti ^ 0.25. 
For the ASCA/ RXTE model s pect rum in the case of a ho- 
mogeneous corona (see Section 3.4), we find the ratio of the 
hard flux, _Fh (without including reflection) to the soft flux, 
Fs, of Fii/Fs ~ 0.50. (The difference between this ratio and 
£h/is is due to scattering of soft photons in the corona.) 
This corresponds to Lh/I/s = 2( _Fh/-Ps) cost (see Appendix 
B). At i = 35° (see Section O), Lh/is « 0.82. Equation 



(B8) yields then / « 0.63. Thus, the geometry with a corona 
covering most of an inner disc is possible in the soft state. 

On the other hand, the corona may be patchy (Haardt, 
Maraschi & Ghisellini 1994; Stern et al. 1995; Poutanen 
& Svensson 1996). We have investigated that possibility 
by adding an additional blackbody disc component to our 
model. The emission of that component corresponds to the 
fraction of the disc photons that are not passing through the 
corona. The remaining fraction, g' , corresponding to the disc 
photons that are incident on the hot plasma, is q' > g (where 
g the covering fraction deflned in Section i.4) because the 
active regions of the patchy corona can be elevated above the 
disc surface. We have found that the addition g' as a free 
parameter improves the flt by Ax^ = —9, i.e., it is statisti- 
cally signiflcant. At the best fit, at which g' = 0.55, the fiux 
ratio is Fh/Fs = 0.53. This then gives, from equation (B8), 
/ = 0.46, relatively close to / in the homogeneous-corona 
case in the paragraph above. 



-6 at Hn = 17tg^. Other disc parameters are now 3.10 Luminosity and accretion rate 



fcTmax = 375Ig eV and (at i = 45° and before the GR cor- 
rection) Mx ~ 1017*^1^40. In contrast to the rin — 6 case, 
no additional constraints on the disc inclination can now 
be obtained from the X-ray data. However, our PN disc 
model does not take into account relativistic corrections to 
the shape of the observed disc-blackbody spectrum. Thus, 
this result should be treated simply as a demonstration that 
the inner radius of a blackbody disc can, in principle, be 
constrained by X-ray spectral fits. 

On the other hand, evolution of power spectra. X-ray 
time lags and the coherence function indicate that on 1996 
May 30 (when the ASCA/RXTE observation fitted above 
was performed) the source was still in transition from the 
hard to the soft state, and the soft state was established only 
later, in particular during the RXTE/OSSF, observations of 
1996 June 17-18 (Cui et al. 1997). This would then support 
our finding of rin > 6 during the ASCA/RXTE observation. 



The total luminosity in the soft state is, using the 
ASCA/RXTE spectrum (see Appendix B), 



is = 47r_D^ 



Fh(l + afl/27r) 



2 cos i 



(9) 



10"* erg cm"^ s"^ 



yield, at i = 35°, L « 4.5 x 10^^ erg 
s~^. This corresponds to ~ 0.03Le at Mx = IOMq. Here, 
Le = 47r^eG'Mxcmp/crT is the Eddington luminosity, where 
/ic ~ 2/(1 + X) is the mean electron molecular weight and 
X ^ 0.7 is the hydrogen mass fraction. 

If we assume rin = 6 and no advection, this corresponds 
to m = Mc^/Le ^ 0.5, or M « 8 X lO" g s~\ On the other 
hand, M can be obtained from fit ting the disc model at 
known Mx and Ls/L, see equation (A7). As expected, this 
yields the same Ad as above. 



4 DISCUSSION 
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4.1 Geometry and scaleheight of the corona 

We consider here constraints on the geometry of the X7 
source from the plasma compactness. For the model spec- 
trum for the ASCA/RXTE data, Lh « 1.9x 10^'' erg s"^ (see 
above). The hard compactness is constrained to 2 £h ^ 7, 
whi ch follows from 5 20 and £h/is~ 0.36 (Section 

3.5). From equation we find the characteristic radius of 
the hot corona is then i?h ~ (H/gRh) (1.2-4) x 10^ cm. At 
the black- hole mass of IOM0 (see Section 3.5), this corre- 
sponds to i?h « (8-30)(Jf/sri?h)i?g- 

We can determine H if we assume hydrostatic equilib- 
rium of the corona. The pressure in the corona has contri- 
butions from radiation, thermal and non-thermal electrons 
and ions, and magnetic field. We find that if the ion temper- 
ature, Ti, equals T^, the dominant contribution to pressure 
in the corona is from radiation. We use then equation (45) in 
SZ94 to estimate H/Rh ~ 10~^ in the case of a homogeneous 
corona {g = 1) at the radius of maximum disc dissipation 
per log R, 17 Rg. This would imply a rather small size of the 
plasma, Rh ^ 0.3g~^Rg, violating our assumption of a ho- 
mogeneous corona, and thus implying g <^ 1. On the other 
hand, Ti 2> Tc is a likely condition in many astrophysical 
situtations, e.g., two-temperature coronae are considered by 
Di Matteo, Blackman & Fabian (1997). In addition, ions can 
be accelerated to non-thermal energies, in addition to accel- 
erated electrons. Then, the corona may be supported by the 
ion pressure, together with the pressure of an equipartition 
magnetic field. 

We can constrain the average energy of ions by requir- 
ing that the rate of Coulomb energy transfer from the ions 
(both thermal and non-thermal) to electrons corresponds a 
luminosity, i'couiomb ^ ^h, which follows from our fits, in 
which Ith <C £h- This, via the condition of hydrostatic equi- 
libirum, constrains in turn the scaleheight of the corona. We 
then obtain, using equations (1) and (3) in Zdziarski (1998), 



H_ 

Rh 



2 X lo^ee^Vth' 



r ic 



(10) 



where Oc = kTc/rricC , and which holds independent of ra- 
dius. In the second equality, we used the numerical values 
of the parameters corresponding to the plasma in Cyg X- 
1. Then, e. g., H / R\^ = 0.3 (corresponding to £th/^h ~ 0.3, 



see Section 



3.5 ) would correspond to 7?h 
0.5 (Section 



(6-20)i?g at the 
. Thus, the hot 



covering fraction of g 
plasma can extend over a substantial part of the disc re- 
gion where most of the gravitational energy is dissipated. 
We note that H/Rh of equation ( |lo| ) corresponds sub-virial 
ions, kTi/mpC^ ~ 0.05/r. 

On the other hand, the corona is not necessarily in the 
hydrostatic equilibrium (e.g. Beloborodov 1999). Some of 
the power released in the corona can be used to accelerate 
it to some velocity, /3c. We can constrain (3 by considering 
the power (or kinetic luminosity), Lkin, of the outflowing 
ions, which equals to the number flux times the bulk-motion 
energy times area. This yields 



(11) 



For r ~ 0.25 and 2 ^ 4 7 (Section (3 ^ 

0.3(I/kiii/ih)^^''. Since no apparent effects of eventual dissi- 
pation of Lkin are observed from Cyg X-1, probably Lkin ^ 



Lh. Then, equation (|l^) allows for a mildly relativistic out- 
flow, e.g., /3 ~ 0.2 at Lkin ~ 0.2Lh and 4 = 7. 

4.2 Coulomb interactions and plasma heating 



As found in Section the ASCA/RXTE data require 
7^1^ 27. The upper limit is due to effects e"'"e~ pair 
production and annihilation, and the lower limit is due to 
energy loss of non-thermal electrons by Coulomb scattering. 

The origin of the lower limit has been discussed by 
Poutanen (1998, see also Zdziarski et al. 1990). For com- 
pleteness, we discuss it in some detail here. The energy loss 
rate due to Compton scattering is (e.g. Rybicki & Lightman 
1979) 

4 ctt , 2 
3 rricC 



'yCompton 



;(7'-l)C/ph, 



(12) 



where Uph is the energy density of photons for which Comp- 
ton scattering is in the Thomson limit. For a homogeneous, 
optically-thin, spherical source of radius TZ and with lumi- 
nosity £, the energy density is roughly 

„2 D 



c 



n 



3 niaC 



ph 



(13) 



(4/3)7^7^3 c 47r (TT 7^' 

where we use I ^ 1^ + (.h for the soft state of Cyg X-1, 
in which most of photons in the observed spectrum are in 
the Thomson limit. On the other hand, the Coulomb energy 
loss rate due to interaction of an electron of 7 ^ 1 and 
7 — 1 ^ Be with a background thermal plasma is given 
by (Gould 1975; Frankel, Mines & Dewar 1979; Coppi & 
Blandford 1992) 

3 c 



'^Coulomb 



271 



rthln(A7i/2) 



(14) 



where A = [27^/(3aoTth)]^/^ and ao = 5.29 x 10"^ cm is the 
Bohr radius. Then, two rates are equal at a critical Lorentz 
factor given by 

37rrth 



7c: 



21 



■InA-f 1. 



(15) 



Note that 7cr depends on IZ only weakly via In A, and thus 
this results is nearly-independent of the source geometry. 
Since 7couiomb oc 7° while 7compton oc 7^, the Coulomb 
energy loss is so fast that it creates a break at 7cr in the 
non-thermal distribution. This then results in a break in 
the distribution of scattered photons, which is not seen in 
our data. Also, the energy lost by the non-thermal electrons 
strongly heats the thermal plasma, increasing its temper- 
ature. These effects are especially pronounced for Finj > 2 
(as is the case in our models), when most of the non-thermal 
energy is provided to accelerated electrons at the low-energy 
end of their distribution. 

Since 7crit oc , the observed lack of signatures of 

the above effect constrains I from below. At the lower limit of 
^ ~ 7, and at rth ~ 0.25, TZ = 10*cm, 7cr ~ 2, i.e.. Coulomb 
interactions just begin to affect the low-energy end of the 
non-thermal distribution. Figure |^ illustrates the effect of 
varying compactness on the spectral shape. 

In most our fitted models, we found Coulomb heating of 
the thermal plasma sufficient to maintain the plasma tem- 
perature required by the observed spectral shape, fcTc ~ 50 
keV (which, we note, is much higher than the Compton tem- 
perature, ~ 1 keV, at which the net energy transfer via 
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Compton scattering between electrons and photons of the 
observed spectrum w ould be null). Thus, Inth/t-h ~ 1 in 
those models (Section 3.5). However, this result depends on 
7min, which has been set so far to 1.3. A higher value of 7min 
would decrease the Coulomb heating and would require some 
additional heating of the thermal plasma component. For 
example, the ASCA/ RXTE data fitted with 7min = 3 gives 
^nth/^h = 0.62, i.e., additional heating is required. However, 
that this fit is worse by Ax^ ~ +9 with respect to the fit 
with 7min = 1.3. Thus, 7min is coustraiued by the data to be 
- 1. 

Some additional heating of the thermal electrons can be 
due to the Coulomb energy transfer from energetic ions, see 
Section The ions can be thermal or non-thermal. Accela- 
ration of ions is, in fact, quite likely, and not in conflict with 
the data. In Section p. 3, we found that the electron accel- 



eration rate, Mnj, is soft, with most of the accelerated elec- 
trons near 7min ~ 1. If acceleration of ions follows a similar 
law, most of accelerated ions will have energies close to the 
low-energy end of the power law, and effects of pion and pair 
production by ions will be negligible. On the other hand, the 
accelerated ions will provide some contribution to Coulomb 
heating of the thermal electrons (and to the presure in the 
plasma), which contribution is allowed by the data as long 
as the average ion energy is ^ 50MeV/r (see Section [4.l[ ). 
We also note that the energetic ions would transfer some en- 
ergy to the non-thermal electrons via Coulomb scattering. 
However, this rate, inversely proportional to the ion mass, 
is found by us to be negligible compared to the energy loss 
rate of non-thermal electrons to the thermal ones [equation 
&]■ 

Another mechanism of heating the thermal electrons 
is the so-called synchrotron boiler (Ghisellini, Guilbert & 
Svensson 1988; Ghisellini, Haardt & Svensson 1998; Svens- 
son 1999). In this process, thermal electrons are heated by 
synchrotron reabsorption of photons emitted by the non- 
thermal electrons. 



4.3 Anisotropy effects 



As stated in Section 3.4, our Comptonization model assumes 
all the distributions (photons and e*) to be homogeneous 
and isotropic. This approximation may not hold in some 
geometries. In particular, in the disc-corona geometry, the 
seed photons enter the corona from below only. Thus, their 
distribution is anisotropic. Then, a soft photon emitted by 
the disc surface is much more likely to scatter on an ener- 
getic electron moving down rather than up. This is because 
the scattering probability is oc 1 — («/c)cosi9, where v is 
the electron speed and i9 is the angle between the directions 
of the electron and the photon. This anisotropy affects in 
different ways spectra from non-thermal and thermal plas- 
mas (and their distributions), and thus we discuss them 
separately. 

The case of a relativistic, non-thermal distribution 
is studied by Ghisellini et al. (1991). Their results show that 
the angle-averaged downward flux, due to a single Compton 
scattering by isotropic electrons with 7 3> 1, is ~ 3 times 
larger than the flux emitted upwards at the angle of 45°. 
This would result in a corresponding strong enhancement of 
the Compton-reflected component, which effect is clearly not 
seen in our data. This enhancement could be much smaller 



if the hot plasma forms an inner hot disc and emission of 
seed photons and reflection are from an outer optically-thick 
disc. 

On the other hand, the steady-state non-thermal elec- 
tron distribution in a disc corona is unlikely to be isotropic. 
This is because electrons directed downward lose energy 
much faster than those directed upward, due to precisely 
the same effect as that discussed above. The main mech- 
anism for their isotropisation is Coulomb interaction since 
higher-order Compton scatterings are negligible in our case 



with Tnth 



10" 



The process of Coulomb interactions of 



non-thermal electrons wi th t he background thermal plasma 
is considered in Section 4.2 above, where we find that the 



observed spectrum constrains compactness to a range where 
the electron energy losses due to this process are negligible 
compared to those due to Compton scattering. Then the 
negligible energy losses, in the case of e*e^ scattering, imply 
no significant changes of the direction of the relativistic elec- 
trons. (The rate of interactions of relativistic electrons with 
protons, which could change the electron direction without 
significant energy loss, is smaller than that with by a fac- 
tor mo/jTip and thus even less signiflcant.) Thus, the steady- 
state distribution of electrons with 7 > 7cr 2, equation 
(p^J will be strongly anisotropic due to the Compton energy 
loss rate depending on the direction. 

When this is the case, an isotropic accelera- 
tion/injection of relativistic electron will lead to isotropic 
emission of singly-scattered photons. This is because the ef- 
fects of anisotropic emission and anisotropic energy loss can- 
cel each other, and the rate of emission in a given direction 
is simply proportional to the injection rate (see Ghisellini et 
al. 1991). 

Note that an anisotropic electron distribution corre- 
sponds to a current in the plasma. This current will be bal- 
anced by an opposite current in the thermal component of 
the electron distribution. Since rth 2> Tnth, the latter current 
will have a negligible effect on the properties of the thermal 
plasma component. 

Anisotropy effects in the case of a thermal distri- 
bution are considered by Haardt (1993), Stern et al. (1995) 
and Poutanen & Svensson (1996). They are due to the same 
physical effect as described above, and thus they suppress 
a part of the emitted spectrum due to the first (and to 
some degree, second) order of Compton scattering. On the 
other hand, photons become closely isotropic after a higher 
number of scatterings. This leads to the so-called anisotropy 
break close to the peak energy of the second-order scatter- 
ing. In contrast to the non-thermal case, we have found that 
the Coulomb relaxation time for the thermal electrons inter- 
acting with themselves is about equal to the thermal plasma 
cooling time (~ 10~* s), and thus no significant anisotropy 
of the thermal electron distribution is expected in our case. 

In order to assess the significance of the thermal 
anisotropy break in the case of Cyg X-1, we have performed 
Monte-Carlo simulations for a slab geometry assuming the 
plasma parameters as found from data fitting. We have 
found the anisotropy break at ~ 15 keV, with the outward 
flux at 7 keV suppressed by a factor of ~ 0.8 with respect to 
the isotropic case. As seen in Figure ^ (the dashed curve), 
this break energy actually coincides with the break occuring 
due to intersections of the power laws due to thermal and 
non-thermal scattering. Then, plasma parameters slightly 
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different from those determined neglecting anisotropy ef- 
fects could still produce the observed spectrum. Thus, the 
present data are still consistent with the presence of a ther- 
mal anisotropy break. 

On the other hand, if the mil dly relativistic outflow is 
present in the corona (see Section |4.l|) , the seed soft radia- 
tion emitted by the disc is more isotropic in the rest frame 
of the outflowing plasma, thus reducing the amplitude of the 
anisotropy break (see Beloborodov 1998, 1999). 



4.4 Disc stability 

We examine here the stability of the optically-thick disc, tak- 
ing into account the fraction, /, of the gravitational energy 
being dissipated in the corona (SZ94). A possible mecha- 
nism for energy transfer from the disc to the corona are mag- 
netic flares (e.g. Galeev, Rosner & Vaiana 1979; Beloborodov 
1999). A part of the disc becomes radiation-pressure domi- 
nate and thus unstable (both thermally and viscously) above 
a critical accretion rate, rhcr. We find that a radiation- 
pressure dominated zone in a PN disc appears first (with 
increasing m) around R = 15.1-Rg. We then extend results 
of SZ94 (who considered the Newtonian potential) to the 
case of the PN potential, and derive nicr in the presence of 
external dissipation as. 



100 



rhcr = O.64(aMx/M0)-'/*[(l - /)(] 



-9/8 



(16) 



where a is the viscosity parameter and is a parameter in 
the vertically-averaged radiative transfer equation, which is 
set by various authors between 0.5 and 2. We follow here 
Shakura & Sunyaev (1973) and set C = 1 (which value is 
in approximate agreement with calculations taking into ac- 
count the vertical disc structure, A. Rozanska, private com- 
munication). We note that relativistic effects, as approxi- 
mated by the PN potential, have a stabilizing effect on the 
disc, i.e., they increase the value of mcr. Specifically, the 
numerical coefficient in equation ( p^ is only 0.37 in the 
Newtonian case (SZ94). 

At Mx = IOM0, / = 0.63, mcr ~ 1.47q~^''*. This 
implies that the d isc in the soft state of Cyg X-1, with rh ~ 
0.5 (Section 3.1C), is stable at any a. Figure 111 shows the 
solutions of the optically thick accretion disc with the corona 
for 4 values of /. 

On the other ha nd, the data do not require rin = 6, 
and a fit in Section p.q gives nn = 17^^''. This solution 
corresponds to a higher m and a lower Mx than in the case 
of Tin = 6. Then, m may possibly exceed mcr and the disc 
instability may be responsible for the large rin.The precise 
value of m in this case depends on the unknown partition of 
the observed luminosity into the fractions from the corona 
above the disc and an inner hot flow at r < rin, possibly 
advective (Narayan & Yi 1995; Abramowicz et al. 1995). 
Therefore, the current data do not allow us to distinguish 
between the possibilities of = 6 and > 6. 



4.5 Alternative models 

One alternative model to the non-thermal Compton model 
presented here is thermal Comptonization, which was pro- 
posed to model spectra of the soft state of Cyg X-1 (Pouta- 
nen et al. 1997; C98; Esin et al. 1998) and other black-hole 
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Figure 11. The solutions of the optically thick disc (Abramowicz 
et al. 1988) with a fraction of accretion power, /, dissipated in the 
corona (SZ94), for Mx = IOMq and a =0.1 at the radius of the 
most significant instability oi R = 15.1-Rg. The curves are labeled 
by their value of /, with / = 0.63 corresponding for Cyg X-1. The 
solid and dashed horizontal lines correspond to the accretion rate 
of Cyg X-1 in the soft state and to the Eddington luminosity at 
an efficiency of r; = 1/16, respectively. 



binaries (Miyamoto et al. 1991; Esin et al. 1997) However, 
as discussed in Section 3.3, this process is unable to account 



for the observed steep power-law tails extending to high en- 
ergies. 

Another possibility for the production of the power-law 
spectra observed in the soft state is that of bulk-motion 
Comptonization (see, e.g., Payne & Blandford 1981; Colpi 
1988; Laurent & Titarchuk 1999). If there is quasi-spherical 
accretion near the black hole, then cold infalling electrons 
could acquire substantial velocities v ^ c, and Comptoniza- 
tion using the large bulk inflow velocity of the electrons (as 
opposed to their assumed smaller thermal motions) could 
give rise to power law spectra like those observed in the soft 
state. This model predicts a sharp cutoff in the spectrum 
at a few hundred keV. The OSSE data, on the other hand, 
clearly show no cutoff until at least 800 keV which strongly 
rules out the bulk Comptonization hypothesis. 

Also, only a fraction, <^ 1, of the disc blackbody pho- 
tons (emitted semi-isotropically with maximum emission per 
log r at r '-^ 20) will be incident on the inner spherical inflow 
(see Figures 1-2 in Laurent & Titarchuk 1999). This appears 
to lead to a normalization of the tail spectrum with respect 
to the disc blackbody much below the observed normaliza- 
tion. Furthermore, the presence of the Compton reflection 
with a covering factor of ~ 0.7 puts strong constraints 

on the geometry of the accretion flow. Such a covering factor 
is natural in a disc-corona geometry, but it appears incom- 
patible with the bulk Comptonization scenario (Figures 1-2 
in Laurent & Titarchuk 1999). 
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4.6 The nature of state transition 

A remarkable feature of state transitions in Cyg X-1 is tliat 
the associated ciiange in tiie bolometric luminosity is very 
small. We have found L « 4.5 x 10^^ erg s~^ in the soft state 



(Section 3.1C), whereas the maximum bolometric luminosity 



observed in the hard state is ~ 3 x lO^'^ erg s~^ (G97; at D = 
2 kpc). If there is no advection, the corresponding change of 
M is small. In the simplest scenario, we expect no advection 
in the soft state provided the optically-thick disc extends to 
the minimum stable orbit. On the other hand, a commonly 
used model for the hard state comprises a hot inner disc 
(Poutanen, Krolik & Ryde 1997; Poutanen & Coppi 1998; 
Poutanen 1998; Esin et al. 1998), in which case advection 
is likely. If the transition from the hard state to the soft 
one is associated with an increase of M, both that increase 
and the advected fraction of the power in the hard state 
have to be very small (Esin et al. 1998). On the other hand, 
advection may dominate over cooling in the hard state, in 
which case the hard-to-soft transition would be associated 
with a decrease of M. 

In either of the above cases, the hard-to-soft state tran- 
sition appears to be associated with a decrease of the min- 
imum radius of the optically-thick disc, Rin- When _Rin is 
large, the inner disc is hot and it receives a major fraction 
(^ 0.9) of the available gravitational energy. The spectrum is 
hard, F ~ 1.6, and the covering angle of the reflector is rela- 
tively small, fl/2n ~ 0.3. At a certain moment, the cold disc 
starts to move inwards, taking in more and more accretion 
power, while the hot phase diminishes. In the final config- 
uration the cold disc extends down to near the last stable 
orbit, while the hot (non-thermal) plasma probably forms 
a corona above that disc, receiving less of accretion power 
(~ 0.6) than in the hard state. This results in a softening of 
the spectrum, to F ~ 2.5, and an increase of the strength of 
the Compton-reflection component, to f2/27r ~ 0.7. 

Another possible model for the hard state has recently 
been proposed by Beloborodov (1999). In that model, the 
X7 source forms a corona above optically-thick disc extend- 
ing close to the minimum stable orbit even in the hard state. 
The hardness of the spectrum and the weakness of Comp- 
ton reflection in that state are then explained by a mildly- 
relativistic bulk motion away from the disc. The hard-to- 
soft state transition would then involve the corona becoming 
static. Possibly, this may be due to the coronal plasma con- 
sisting of e"'"e~ pairs in the hard state (which is allowed by 
the form of hard-state spectrum, G97; Poutanen & Coppi 
1998), and of an ordinary plasma in the soft state. The 
data studied here appear not to allow distinguishing between 
these possibilities. 



5 CONCLUSIONS 

In this work, we have obtained a self-consistent descrip- 
tion of the soft state of Cyg X-1 (observed in 1996 by 
ASCA, RXTE and OSSE) in terms of a (predominantly) 
non-thermal corona above the surface of an optically-thick 
accretion disc. Emission of the non-thermal corona accounts 
for both the high-energy power-law tail observed up to sev- 
eral hundred keV and the soft excess, resolving two problems 
of the thermal-plasma model of C98. 



The non-thermal corona receives most of the supplied 
power in the form of electrons accelerated (or injected) to 
relativistic energies. The acceleration rate is rather soft, with 
the power-law index of Finj ~ 2.5-3 above 7min ~ 1. The 
electrons lose their initial energy in Compton and Coulomb 
processes, and flnally form a thermal component (at kTc ~ 
50 keV) of their distribution at low energies. The resulting 
X7 spectrum consists of blackbody photons emitted by the 
disc (at low energies) and a component due to Compton 
upscattering of the disc photons by both thermal and non- 
thermal electrons in the corona. 

Production of e'^e~ pairs is negligible in the corona, 
which follows from the observed lack of an annihilation fea- 
ture in the spectrum from OSSE. This constrains from above 
the compactness (a ratio of the luminosity to the size), £, of 
the source. On the other hand. Coulomb scattering is an 
important process in heating the thermal electrons by the 
non-thermal ones. Still, the Coulomb energy loss by non- 
thermal electrons is constrained by the observed spectrum to 
be slower than their Compton energy loss, which constrains 
£ from below. The resulting constraints are 7 ^ £ ^ 30. 

We also flnd that a modest heating of the thermal elec- 
trons (in addition to their Coulomb heating by the non- 
thermal electrons) is required by the observed spectrum. 
This heating can be provided by Coulomb interactions with 
energetic ions in the corona, which then also provide most 
of the pressure support of the corona. 

The coronal geometry of the source is supported by the 
observed spectral component from Compton reflection from 
the disc surface and the Fe broad Kq fluorescence line. Both 
components are relativistically smeared, supporting reflec- 
tion from an inner part of the disc. The strength of the ob- 
served Compton reflection spectrum is fully consistent with 
reflection from a slab, after taking into account Compton 
scattering of reflected photons in the coronal plasma. The 
observed softness of the spectrum is consistent with energy 
balance in the disc-corona system with / ~ 0.6 of the avail- 
able gravitational energy dissipated in the corona (and the 
remainder within the disc). From the observed luminosity 
and fitted compactness, we find the corona covers a sub- 
stantial fraction of the disc provided most of its pressure 
support is provided from energetic ions. Still, the ion en- 
ergy in random velocities is constrained to be much below 
the virial energy (by the constraint on the Coulomb heating 
of thermal electrons, see above). Then, the corona can be 
static. On the other hand, we find that a mildly relativistic 
outfiow is also consistent with the observations. 

Fits of our PN disc model together with the mass 
function and constraints on the mass ratio imply Mx — 
(10 ± 1)M0 and the inclination of i ~ 30°-40°, at the 
most likely distance of 2 kpc and the colour correction of 
/col = 1.7. The bolometric luminosity is ~ 4.5 x lO^'^ erg 
s~^ ~ 0.03I/E, which corresponds, in the absence of advec- 
tion, to m ~ 0.5. At that accretion rate, the disc is fully sta- 
ble all the way down to the minimum stable orbit against the 
instabilities associated with the dominance of radiation pres- 
sure. On the other hand, solutions with advection, a higher 
m and the disc terminating at rin ^ 6 are also allowed by 
the data. 
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APPENDIX A: SPECTRUM FROM A 
PSEUDO-NEWTONIAN DISC WITH 
ZERO-STRESS INNER BOUNDARY 

Let us consider a geometrically-thin, optically-thick, accre- 
tion disc around a Schwarzschild black hole with a mass, 
M. We assume the disk locally emits a blackbody spectrum, 
which fraction, p^c is then partly scattered in a corona. The 
unscattered flux observed from the disc is given by 



2-K cos i 
fi n2 



B^{T)RAR, 



(Al) 



where Bi, is the blackbody spectrum. Other symbols used in 
this Appendix are defined in Sections |l| and 3.1 except that 
we now drop subscripts in symbols for the colour tempera- 
ture and mass. We allow here for truncation of the disc at 
> R^^ = &GM/c^. 

In order to compute the observed spectrum, we need to 
find the radial colour temperature distribution, T{R). The 
rate of energy generation per unit area of one side of the 
disc, Q, is given by 



Lu{Rn-is)Rn 



LOR^ 



(A2) 



where uj is the rotation velocity. The colour temperature is 
given by. 



qQ{R) = <t\^ 

/col 



(A3) 



where a is the Stefan-Boltzmann constant and q is the ratio 



of the disc flux to Q{R), see equation (B6) below. When 
a fraction (/) of Q{R) is dissipated in a corona above the 
disc surface (SZ94), g < 1, and the disc emission is due 
to both internal dissipation of (1 — f)Q{R) and irradiation 
by the corona. We adopt the pseudo-Newtonian potential 
(Paczyhski & Wiita 1980), 



GM 
R - 27?, 



(A4) 



The rotation velocity is then 



oj{R) 



Then we find from equations (AS)-(A_2) 



T{r) = To 



r - 2/3 



r(r - 2)3 
where r = R/Rg and 
3Mqc^' 



33/2 



21/2 ^3/2 



1/4 



To — fa 



1/4 



(A5) 



(A6) 



(A7) 



■Tor 



-3/4 



The maximum T is 



(A8) 



At r > 6, T(r) 

r CqTo, rin < I max , 

\r(nn), otherwise, 

where co ~ 0.1067 and rmax ~ 9.505. 

The (unscattered) flux emitted by the disc per u nit p ho- 
ton energy, E, can be then obtained from equation (Al), 

rdr 



Fe = KE-' 



where 



K ^ {1 



exp[E/kT{r)] - 1' 



An / GM \ ^ cos i 



'/l3c2 Vc2Dy 



J-4 
J col 



(A9) 



(AlO) 



and h is the Planck constant. This model, under the name 
diskpn, is a part of the xspec packet. Its parameters are 
Tmax, rin and the normalization, K. 

F rom these parameters, M can be obtained from equa- 
tion ([A10|), and then M can be obtained using equations 



For a constant K, M oc Df^^^ cos 



-1/2, 



For con- 



stant K and Tmax, M oc f^^i cos 
is given by 



Ls = 



FeAE, 



The disc luminosity 



(All) 



which, for 



- 6, equals 

1 . 2 

-Psc)«7 — Mc 



O'^O./col'^ 



(A12) 



Finally, we point out that the above model differs signif- 
icantly from a commonly-used disc blackbody model, diskbb 
(Mitsuda et al. 1984) in xspec, in which the torque-free 
boundary condition is neglected and the asymptotic, r S> 6, 
distribution is assumed, i.e., T{r) = Tor~^^'^. This is a suffi- 
ciently accurate model for nn S> 6. However, when the disc 
extends close to the minimum stable orbit, the discrepancy 
between this T{r) and the T{r) of equation (AC) becomes 



very significant, as shown on Figure Al. Still, a diskbb spec- 



trum can be made relatively accurate with some rescaling of 
-Rin and Tmax- We find in the case of rin = 6 that a PN spec- 
trum at some Rin and Tmax is close to the diskbb spectrum 
computed for -Rf^""'"' « 2.737?in and T^l=f ^ 1.04Tmax. 

A more serious discrepancy regards the accretion effi- 
ciency. The efficiency (77 = L/ Mc?) of a PN disc extending 
to the minimum stable orbit is 77PN = 1/16, very close to the 
exact value in the Schwarzschild metric, r;s = 1 — (8/9)^''^ ~ 
0.0572 (e.g. Kato, Fukue & Mineshige 1998), whereas the 
diskbb model yields r^diskbb = 1/4. Therefore, M for discs 
extending to the minimum stable orbit computed using the 
diskbb model and assuming t^s (as sometimes done) is 4.37 
times less then the actual M . 
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Figure Al. The radial te mpe rature distribution in a pseudo- 
Newtonian disc, equation ( [Aq ) (lower curve) compared to the 
asymptotic distribution, oc r"''/'', used by the diskbb model (up- 
per curve). 



APPENDIX B: ENERGY BALANCE IN 
DISC-CORONA SYSTEMS 

We express here energy balance between a corona and an 
underlying disc (see Figure^ in terms of the observed soft, 
hard and reflected fluxes, F^, _Fh and F^, respectively. This 
allows a direct comparison with results of data fitting, in 
which we can distinguish between the flux in unscattered 
disc photons (Fs), the flux in scattered photons (-Fk), and 
the flux in the unscattered Compton-refiection and Fe Ka 
photons (Fr). 

We first note that Fs and fh are related in different 
manners to the corresponding luminosities. The disc is opti- 
cally thick and thus it radiates with a semi-isotropic specific 
intensity, which implies Lb = {2-k / cos i)D^ Fs. In general, 
the plasma emission can be anisotropic as well. If the emis- 
sion radiated by the corona in directions away from the disc 
is isotropic (which we assume in this work), Lh = AtxD^F\i. 

We consider then a corona dissipating a fraction, /, of 
the total available power, with the underlying disc dissipat- 
ing the remaining (1 — /) fraction (SZ94). In general, the 
corona can be either homogeneous or patchy. In the latter 
case, Z/B includes also emission of the part of the disc not cov- 
ered by the corona. We will denote the powers dissipated in 
the corona and in the disc by Pc (= fL) and Pd [= (1 ^ /)^] , 
respectively. From energy conservation. 



L = Pc + Pd = + + 



(Bl) 



We will allow here for an anisotropic corona emission, in 
which the power emitted down equals d times the power 
emitted up. Also, we will include the power of Compton- 
refiected photons scattered in the corona in the total emis- 
sion of the corona, (1 -|-d)I/h. Thus, the total emission of the 



corona includes its own dissipation and the power in photons 
emitted by the disc and then scattered, 

(l-|-d)Lh = Pc+Psc(Pd + dih), (B2) 

where psc is the average probability of scattering of a photon 
emitted by the disc. This can be solved for 

Pc + PscPd 



d(l 



(B3) 



The luminosity in unscattered soft photons emitted by the 
disc includes parts due to the intrinsic disc emission and due 
to the reemission of the incident coronal emission. 



Ls = (l-Psc)[Pd + d(l-a)Lh] 



(B4) 



(where a is the integrated albedo of the disc). The above 
equations can be solved for 



Le = 



(I - Pec) {[1 + d(l - apec)] Pd + d(l - a)Pc} 



(B5) 



l + d(l-psc) 

The ratio of the soft fiux at the disc surface to the total 
dissipation rate per unit area is 

Lb 



(1 - Psc)L 



(B6) 



(see Appendix A). Finally, the luminosity in unscattered 
Compton-refiected photons is 



Lr = (1 - Psc)dfflih. 
The solution for / is 

1 + d(l - apsc) - Psc(l 



/ = 



^Ls/Lh 



(B7) 



(B8) 



1 -I- da(l - Pbc) + Ls/Lti 

We note that our results are consistent with those of 
Poutanen & Svensson (1996). However, their condition of 
energy balance is expressed in terms of the total soft emis- 
sion of the disc (denoted by Zs), which includes the part scat- 
tered in the corona, and the total emission of the corona (Zc), 
which includes the (unobserved) coronal emission directed 
towards the disc surface. For comparison with observations, 
those quantities need be transformed to the luminosities cor- 
responding to observed fiuxes; e.g., Zc corresponds to (1-l-d) 
times the actual luminosity of the corona. 

This paper has been produced using the Royal Astronomical 
Society/Blackwell Science MJi?^ style file. 
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